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Abstract 

We consider the electromagnetic (EM) perturbative effects produced by the high-frequency 
gravitational waves (HFGWs) in the GHz band in a special EM resonance system, which 
consists of fractal membranes, a Gaussian beam (GB) passing through a static magnetic field. 
It is predicted, under the synchroresonance condition, coherence modulation of the HFGWs to 
the preexisting transverse components of the GB produces the transverse perturbative photon 
flux (PPF), which has three novel and important properties: (l)The PPF has maximum at a 
longitudinal symmetrical surface of the GB where the transverse background photon flux (BPF) 
vanishes; (2) the resonant effect will be high sensitive to the propagating directions of the 
HFGWs; (3) the PPF reflected or transmitted by the fractal membrane exhibits a very small 
decay compared with very large decay of the much stronger BPF. Such properties might provide 
a new way to distinguish and display the perturbative effects produced by the HFGWs. We also 
discuss the high-frequency asymptotic behavior of the relic GWs in the microwave band and the 
positive definite issues of their energy-momentum pseudo-tensor . 
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I. INTRODUCTION 



Unlike usual celestial gravitational waves (GWs) having low frequencies, which are often 
a small fraction of a Hz, the relic GWs in the microwave band (~ 10 8 — 10 11 Hz), pre- 
dicted by the quintessential inflationary models (QIM)_l[ [2], 3}, the pre-big bang scenario 



aodels (QIM) 1 
rios fl fl0,3, 



form high-frequency random sig- 



(PBBS) and some string cosmology scenarios 
nals, their root-mean-square (rms) values of the dimensionless amplitudes might reach up 
to ~ 10~ 30 — 10~ 33 / y/llz, and because of their weakness and very high-frequency proper- 
ties, they are quite different from the low-frequency GWs. The thermal motion of plasma 
of stars, the interaction of the EM waves with interstellar plasma and magnetic fields, the 

rf 

evaporation of primordial black ho 



energy laboratory schemes 
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es |a | , even ultra-high-intensity lasers [9J and other high- 



ly are possible means to generate the HFGWs in the 



GHz band and higher frequencies. Interaction of the HFGWs with the EM fields and the 



EM detection o 
aut 



nors 



36 



37 
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: the HFGWs have been theoretical 
15, 3, 17, 3, 3, 20, 21, 22, 23, 24, 



y and experimentally studied by many 



25j, 



26 



27 



28 



29 



30, 



31 



32, 



33 



34 



35, 



.These works include the gravitation-EM conversion in the static EM fields(e.g., 
the Gertsenshtein effect and its inverse effect), the cavity classical- and cavity quantum- elec- 
trodynamical response to the HFGWs, resonant photon-graviton conversion, Barry's phase 
in the EM detection of the HFGWs, resonant interaction of the HFGWs with the EM wave 
beams, the rotation of the polarization vector of EM wave caused by the HFGWs in the 
toroidal waveguide, the difference frequency resonant response of coupled spherical cavities, 
etc. 

Although the relic GWs have not yet been detected, we can be reasonably sure that 
the Earth is bathed in the sea of these relic GWs. Since 1978 such relic and primordial 
background GWs have been of ever increasing scientific interest as many researches have 



shown 



39, 



40 



41 



42j. 



Based on high-dimensional (bulk) spacetime theories, it has also been theoretically shown 



43|, |44j|that all familiar matter fields are constrained to live on our brane world, while 



gravity is free to propagate in the extra dimensions, and the HFGWs (i.e., high-energy 
gravitons) would be more capable of carrying energy from our 3-brane world than lower- 
frequency GWs. It is noted that propagation of the HFGWs may be a unique and effective 
way for exchanging energy and information between two adjacent parallel brane worlds 
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45 



46 



47 



48| . Moreover, if the pre-big bang scenario is correct, then the relic GWs would 



be an almost unique window from which one can look back into the universe before the big 
bane; , l^] . Although these theories and scenarios may be controversial and whether or 
not they have included a fatal flaw remains to be determined. The successful detection of 
the high-frequency relic gravitational waves (HFRGWs) will certainly shed light on many 
of these theories. 

In this paper we shall discuss some ideas and theoretical basis for selection and detection 



of the HFGWs with thepredicted typica 
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-30 



10 



-33 



'Hz 



epreaictf 
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parameters v g ~ 5 x 10 9 Hz (5GHz) and h r 



40]. This paper includes following seven parts: (1) 



Introduction; (2)Asymptotic behavior of the relic GWs in the high-frequency (the microwave 
band) region and positive definite issues of the energy-momentum pseudo-tensor of the 
HFRGWs; (3) The EM resonant system to the HFRGWs, i.e., the coupling system of the 
fractal membranes and a Gaussian beam (GB) passing through a static magnetic field; 
(4)The EM resonant response to the HFRGWs and some numerical estimations; (5) Selection 
and detection of the PPFs; (6) A very brief review to the noise issues; (7) Concluding 
remarks. 



II. THE HIGH-FREQUENCY RELIC GRAVITATIONAL WAVES IN THE GHZ 
BAND 

1. High-frequency asymptotic behavior of the relic GW in the microwave band. 

It is well known that each polarization component hij(rj, x)oi the relic GW can be written 

as A, Em 

hij = exp(ik ■ x)eij, (1) 

(X 

The time dependent of h is determined by the /x(r/) satisfying the equation 

ji+ (k 2 -a/a)n = 0, (2) 
where a = |^§, a = a{rf) is the cosmolo gy s cale factor, rj is the conformal time. In fact, the 



Eq. (j2J) has different exact solutions 5jJ, l52j in the different evolution stages of the Universe, 
their analytic forms are often very complicated. Fortunately, for HFRGWs in the GHz band 
(i.e., the relic gravitons of large momentum), we have k 2 ^> \a/a\ in Eq.(|2]), i.e., term a/a 



can be neglected, then the solution forms can be greatly simplified. In this case Eq. (J2]) has 
the usual periodic solution 

fi(rj) = Ai(k) exp(—ikrj) + A 2 (k) exp(ikrj). (3) 

By using Eqs. ([1]) and ([3]), we have 

h = A 1 (k)/a(rj) exp[i(k ■ x — kr/)} + A 2 (k)/a(rj) exp[i(k ■ x + krj)], (4) 

Consequentially, the HFRGWs can be seen as the superposition of all "monochromatic 
components", Eq. PJ. 

2. The positive definite issues of the energy-momentum pseudo-tensor of the HFRGWs. 

If the relic GWs do exist and have an observable effect, they should have reasonable 
expressions for their energy-momentum pseudo-tensor (EMPT). In particular, the energy 
density of the relic GWs should be positive definite, and the momentum density components 
should have reasonable physical behavior. Although the energy spectrum of the relic GWs 
and their imprint on the cosmic microwave background have been much discussed, there is 
little research into the complete forms of the EMPT of the relic GWs 3), 5^], this research 
may provide a theoretic basis for existence of relic GWs and their detection. Unlike previous 
works, our attention will be only focused on the EMPT of the HFRGWs in the GHz band, 
especially the positive definite property of the energy density of them. 

The relic GWs are small corrections to the background metric tensor, the spacetime 
background is the de Sitter spacetime, and the metric takes the form 
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40] 



ds 2 = a 2 (r])[—dr] 2 + (8^ + h i j)dx l dx :j ] = g^dx^dx" , (5) 

Qtiv = 9fiu + a 2 h^ u , (6) 

where 

gp,v = (-a 2 , a 2 , a 2 , a 2 ), (7) 



Sij is Kroeneker symbol. Because the sea of the HFRGWs can be seen as the superposition 
of all "monochromatic components" of the Fourier expansion, each "monochromatic com- 
ponent", Eq. (j3J), contains every possible propagating direction. In this case, we consider a 

4 



single "monochromatic wave" propagating along the z-axis in Cartesian coordinates without 
loss of generality. From Eqs. (jSJ)- ©, then the metric has following form in Cartesian 
coordinates 



( -a 2 

a 2 (l + /i f 
a 2 h® 
\ 



a 2 (l - hr 



\ 




a 2 / 



(8) 



From Eq. ([H]), we have 



#00 = -a 2 , 9u = a 2 (l+/i©), g 2 2 = a 2 (l-h e ), g 33 =a 2 , g 12 =g 21 =a 2 h i 



(9) 



and 



g = det(g^ v ) = a (h e + - 1) 



(10) 



Expressions of the Einstein EMPT are 541 ] 



where 



TJVCT 



9 » - \QttG 



--9,x[-g(g u Y 7 - 9^)1 



(ii) 



(12) 



is the super-potential. Since h 2 terms have been taken into account in the determinant of 
metric. Eq. fllOp . and the EMPT for gravitational field concerns quadratic terms of h, the hi 
and h 2 ^ terms in the inverse of metric should also be taken into account. By using Eqs. 
fl8|).fl9l. f|T0|) and g^ a g aiJ = 5^, neglecting third- and higher-orders infinitely small quantities, 
we obtain non-vanishing components of g^ u and H^ a in empty space as follows 



9°° = ~a\ 



g u = a - 2 (l-/i e + h% + hi), g 22 = a- 2 (l+/i e + h 2 @ + h 2 



(13) 



g 33 =a 2 , g 12 =g 21 = — a 2 h li 



tt03 

^o 



--9oo (-gg 00 g 33 ) 3 = -^ka 2 (hl + hl 



(14) 



5 



Hi 1 = -H[ 



10 



^9n (g^g^+^g^g^ 



Aaa — ah® — 2aa(h (B + h#) — a (h^h® + h^h®), 



02 



=9n (^ 2 V ) )0 + ^ 12 (^ 2 V°), 



-H 20 = 

a 2 {h^h(Q — h^h^) — a 2 h^, 



Hi 3 = -HI 



51 ' 5ii (-99 1 V 3 ) 3 + -}=m (-99 12 9 33 ) 



-g 



-g 



-ika (h® + h® + h^) 



^23 H 32 



=9n (-99 21 9 33 ) 



--9i2 {-gg 22 g 33 ) 3 = -ika 2 h, 



H 01 = -H 



Hf 



10 



' 9n(99 U 9 00 ). + -^9i2 (99 12 9 °) 



'-9 



-9 



a 2 {h§h^ — — h§ 



-Hf 



^ (99 2 y°), + -^922(99 22 9 °), 



Aaa + a — 2aa(h (B + h\) — a (h^h® + h^h®), 



H\ 3 = -Hf = — £ 21 (-gg n g 33 ) 3 + -^g 2 2 (-gg 12 g 33 ) 3 = -i*a 2 /i { 



-Hf 



-^ 9 2i (-gg 21 g 33 ) 3 + -L^ (-<7<7 2 V 3 ) 

ika 2 (h e -h%- h%), 



rr03 _ rr30 
H 3 - ~ H 3 



'-9 



533 (gg 33 g 00 ) 



Aaa — 2aa{h (B + h^) — 2a (h^h® + h®h$ 



From Eqs. P|)- f|T0|) . f|T3l) and f|T4l) . we obtain the energy density and energy flux 
follows, respectively 

r 4 4 k 2 
j.0 _ C ttOo _ __2_(h 2 _L h 2 \ 

~ 16ttGV=5 0,(7 ~ 4vrGa2 lft ® + ^' 



ct\ = cH 1 ^ = ct 2 = cH { 



0,(7 _ 



(25) 




Tj3(J 



0,cr — 



47rGa 3 



ike 5 



[a{h% + h%) + a(/i e /i e + h^h^)]. 



(26) 



For the "monochromatic components" of propagating along the x- and y-axises, we have the 
same expression for the energy density, Eq. ( |24|) . but t$ = t% = and = ijj = 0, respectively 
The energy flux density also has the same form, Eq. ( 1261) . Thus, the energy density of the 
HFRGWs is positive definite, and the energy flux densities have reasonable physical behav- 
ior in the conformal time coordinates. If we integrate the EMPT for all "the monochro- 
matic components" of the HFRGWs, then we can find that the EMPT is homogeneous and 
isotropic. Riazuelo and Uzan [3!]obtained an expression for the EMPT of the relic GWs in 
the momentum space, the average values of such expressions for the EMPT have reason- 
able physical behavior. However, Eq. (jlj) shows that the stochastic relic GWs background 
contains every possible propagating direction, and because of stochastic fluctuation of the 
amplitudes of the HFRGWs over their bandwidth, detection of the HFRGWs will be more 
difficult than that of the monochromatic plane GWs. In this case, can the HFRGWs be 
selected and measured? In particular, if two HFRGWs have the same amplitude and fre- 
quency, but propagate along the exactly opposite directions (standing wave), will their effect 
be canceled and nullified? We shall show that in our EM system the EM perturbation pro- 
duced by the HFRGWs, which propagate along the positive and negative directions of the 
symmetrical axis (the z-axis ) of the GB, will be non-symmetric and the physical effect 
generated by the HFRGWs propagating along other directions will be also quite different, 
even if they satisfy the resonant condition (u e = u g ), and only the HFRGW component 
propagating along the positive direction of the symmetrical z— axis of the GB can generate 
an optimal resonant response. Thus our EM system design will be very sensitive to the 
propagating directions as well as the frequencies of the HFRGWs, and it may provide a 
HFRGW map of the celestial sphere (similar to the map of the relic microwave background 
provided by the Wilkinson Microwave Anisotropic Probe or WMAP). 
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III. THE ELECTROMAGNETIC RESONANT SYSTEM: THE COUPLING SYS- 
TEM OF THE FRACTAL MEMBRANES AND THE GAUSSIAN BEAM (GB) 
PASSING THOUGH A STATIC MAGNETIC FIELD 



Our EM system consists of the GB of a fundamental frequency mode 55j operating in 
;he GHz band immersed in a static magnetic field, with a new-type of fractal membranes 



56, 



57 



581 ] to focus PPF signal along the detection axis. In order to consider resonant 



response to the HFRGWs in the laboratory frame of reference, all parameters of the 
EM system should be values in the frame of reference. The general form of the GB of a 
fundamental frequency mode [55] is 



(27) 



where r 2 = x 2 + y 2 , k e = 2ix/\ e , f = 7rW 2 /X e , W = W [l + (z/f) 2 ] 1 / 2 , R = z + f 2 /z, ^ is 
the amplitude of electric (or magnetic) field of the GB, Wo is the minimum spot radius, R is 
the curvature radius of the wave fronts of the GB at z, u e is the angular frequency, A e is the 
EM wavelength, the z-axis is the symmetrical axis of the GB, and S is an arbitrary phase 
factor. For the resonant response to a HFRGW, S is the phase difference between the GB 
and the resonant component of the HFRGW. Using a new approach, different from Refs. 
[21, 22], we choose the GB with the double transverse polarized electric modes (DTEM), 
and utilize the coupling effect between the fractal membrane in the GHz band and the GB 
passing through a static magnetic field. Indeed, the GBs with the DTEM exhibit more 
realizable modes, they have been extensi vely discussed and applied in the closed resonant 
cavities, open resonators and free space [55j, |59|, [60|, [6JJ, including the standing- wave-type 



and traveling-wave-type GBs. Moreover, a very important property of the EM system is 



that the PP 



small decay [56 



signal) reflected or transmitted by the fractal membranes exhibits a very 



57 



in transit to the detectors (high-sensitivity microwave photon flux 
receivers) compared with the very large decay (typical Gaussian decay rate) of the much 
stronger BPF. This property provides a possibility to distinguish them in some suitable 
regions. 

If the static magnetic field pointing along the y-axis is localized in the region —l\ < z < I2, 

+ ^ = and 

ox cry 



setting E, 



(0) 



ip = ipx and using divergenceless condition V • E 



B(°) = -i-V x E<°) (we use MKS units), then we have 

Ei 0) = 1> = V,, Ef) =^y = -j ^dy = 2x(^- 2 - i^) J ^ x dy, E® = 0, (28) 
6(0) _ h(o)__±^h m _ i M x _dil) y 



and 



(2 < — /1 and 2 > Z 



2 



where the superscript denotes the background EM fields, the notations ~ and A 
stand the time- dependent and static EM fields, respectively. For the high-frequency EM 
power flux (or in quantum language: photon flux), only non- vanishing average values of 
this with respect to time have an observable effect. From Eqs. fl27|) . fl28|) . and fl29|) . one finds 



(0) 



1(1(^)5(0)) - 1(^)5(0))) 



ZlWe /i A*0 



:M<#>i+<<t>iH o, -<><-^ (33) 



where 7zo; e is the energy of single photon, ni -*, n y ^ and represent the average values of 
the x-, y— and z- components of the BPF densities, in units of photons per second per square 
meter, propagating along the x-, y- and z-axes, respectively, the angular brackets denote the 
average over time, /i , fy^ and /J are the functions of if>o, Wo, ou e , r and z. Because of 
the non- vanishing rix and nf\ the GB will be asymptotically spread as \z\ increases (i.e., 
the irradiance surface of the GB spreads out in the + z and - z directions). 
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IV. THE EM RESONANT RESPONSE TO THE HFRGWS 



For the EM resonant response in the laboratory frame of reference, we should use the 
intervals of laboratory time (i.e., cdt = a{r])drj) and laboratory frequency of the HFRGWs. 
In this case, Eq.(pE]) can be written as 

h(x,t) = A(k g )/a(t) exp[i(k 5 .x— u> g t)] + B(k g )/a(t) exp[i(k g .x.+u> g t), (34) 

where A/ a and B/a are the stochastic values of the amplitudes of the HFRGWs in the 
laboratory frame of reference, k g and u g are the corresponding wave vector and angular 
frequency in the frame of reference. Eq. can be seen as the approximate form of 

each "monochromatic polarization component" of the HFRGWs in the GHz band. In our 
EM system, since only "monochromatic component" of the HFGW propagating along the 
positive direction of the symmetrical axis (the z-axis) of the GB generates an optimal 
resonant response (see Section V), our attention will be focused to a circular polarized 
"monochromatic component" of the HFRGW in the z-direction, i.e., 

= h u = ~ h 22 = A® exp [i (k g z - u g t)] , 

(35) 

h® = hyi = h 2 i = iA® exp [i (k g z - uo g t)} , 

where A^^A® rs A(k g )/a(t) [see Eq.( )34|) ] .Using the electrodynamical equations in the 
curved spacetime 

' <} ^gg m g vP F aP )=^, (36) 



'—g dx v 

V a F^ + V u F afl + V^F ua = 0, (37) 



we can describe the EM perturbation produced by the HFRGWs in the EM system, where 
is the EM field tensor, and F^ = F^ + Ffi* , Fffl and PjiJ represent the background 
and first-order perturbative EM fields respectively in the presence of the HFRGWs. 
indicates the four- dimensional electric current density. For the EM response in vacuum, 
J M = in Eq. (I3"6"|) . Because of the weak field property of the HFRGWs, the perturbation 
methods will still be valid. Using Eqs.(l9l).( fT0i) and (TI3"]) . then Eqs.( l36l) . (1371) can be reduced 
to 
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V. + F«) + V, (f$ + F«) + V, (f£> + *£>) = 0, (39) 

Unlike plane monochromatic GWs, the amplitudes of the relic GW in Eq. (1341) are not 
constant, in this case solving Eqs. (l38|) and (|39|) will often be difficult. In our case, fortunately, 
since this is the EM response in the GHz band, and considering Eq. fl35|) . the following 
equivalent relations would be valid provided u g ^> a/ a, i.e., 

d 

— -> Tiw g , V -> *k 9 . (40) 

In this case the process of solving Eqs.(l38l) and ( |39l) can be greatly simplified without 
excluding their essential physical features. 

Introducing Eqs.flUhCG^CGID and (J35D into Eqs.flMD and (02]), considering |/i e | , \h®\ < 1, 
using the equivalent relations, Eq. (1401) and neglecting high-order infinitely small quantities, 
the first-order perturbative EM fields generated by the direct interaction of the z-component 
of a certain "monochromatic wave," Eq. fl35|) . with the static magnetic field B y ^ can be 
given by Q Q, 28] 



Ex = TjA^By^ k g c(z + h) exp[i(k g z - uj g t)} + ^^Bf^ cexp[i(k g z + u g t)], 

= ^A (B Bf ) k g (z + h) exp[i(k g z - uj g t)] - jA^By^ exp[i(k g z + oJgt)], 
Ey = -^A®By°^k g c(z + h) exp[i(k g z - uj g t)] + jA (S) B y °^cexp[i(k g z + uJ g t)] : 
BP = \A^B i y )) kg(z + h) exp[i(k g z - uj g t)} + {A®B y 0) exp[i(k g z + u g t)], 



(41) 



where A e ,A® m A(k g )/a(t), [see, Eq. ()34l)].— U < z < l 2 . Equation (T4TD shows that the 
first-order perturbative EM fields have a space accumulation effect (oc z) in the interacting 
region: this is because the GWs (gravitons) and EM waves (photons) have the same 
propagating velocity, so that the two waves can generate an optimum coherent effect in the 

l^Ejj ]. Such results and the calculation by Feynman perturbation 



propagating direction 

techniques in Ref.fl^ are self-consistent. In our EM system, we shall neglect the EM 
perturbation solution which describes the EM perturbation propagating along the negative 

11 



direction of the z=axis since it can not satisfy the boundary condition F t 



Obviously, this is typical inverse Gertsenshtein effect 



(i) 



0. 



13|. From Eqs. ([28D , ([29j) , pjj) and 



(|4T|) . the total EM field tensors in the presence of the HFRGW can be written as 



,(o) 



',(1) 



( 

-\{E^ + F 

V o 



y , 



B 



(0) 



(l)^ 



l/p(0) , pi 



(0) 



-5 



bp + m ] 



r(o) , 5(0) , 5(i) 

Dy + Dy + -Dy 







/ 



(42) 



In our exemplar EM system we have chosen the GB power of P = 10W and the static 
magnetic field of = 3T, then corresponding magnetic field amplitude of the GB is only 
£?(°) 10 _5 T, so the ratio of B^ and the background static magnetic field B^ is roughly 
BW/Bjp ~ 1(T 5 . In this case we have neglected the perturbation EM fields produced by 
the directed interaction of the HFRGW with the GB. 

Using the generic expression of the energy-momentum tensor of the EM fields in GW 
fields 

= —(—F£F va + -g^F^F^), (43) 
Ho 4 



we can calculate the perturbation to the energy-momentum of the EM fields in the GW 
fields. Because of the weak field property of the HFRGWs, the energy-momentum tensor 
can also be decomposed into 



(0) (1) (2) 



(44) 



(0) (1) (2) 

where T^ v is the energy-momentum tensor of the background EM fields, T^and T^ v are 

(0) 

first- and second- order perturbations to T^ u in the presence of the HFRGW. From Eqs. 

(0) (1) ' (2) 

(HI and (El, T^T^ and T^ u can be written as 



(i) 

rpflU 



rpilu _ J_[__pM(0) jpua(0) ^i(O) pa(3(0)^ 

Ho 4 a/3 - 

L [-(F£ (0) F va{1) + F^ l) F ua{f)) ) + -^ u (F (l ^F am + F^F aPW ) (46) 



in) 
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(2) 1 11 

Eqs. (pi). (HB|> . and fl4Tj) show that the first-order perturbation of the EM fields 

(i) 

tensor contains only the first-order term of the metric h, thus T^ u is proportional to the 

(2) 

first-order terms of h, while T^is proportional to the second-order terms of h. Because 
the expected amplitude of HFRGWs in the GHz band would be only h~lCT 28 -iCr 33 /V-£^ 



40 



(0) (1) (2) 

52[,then for nonvanishing T Ml/ , T MW and T Ml/ , we have 



(0) (1) (2) 



(1) (0) (2) 

In this case, for the effect of the HFRGW, we are interested in T^ v but not in and T^ v . 
Considering the transverse and traceless (TT) gauge condition (h 11 = —h 22 = h®, h 12 = 
h 21 = h^,h\ = 0, /i 01 = h 02 = h 03 = h 13 = h 23 = h 33 = 0), all novanishing components of 
the first-order perturbation to T^ v generated by a "monochromatic component" propagating 
along the z-axis of the HFRGW can be written as 

(i) i i 

/i 4 af3 

(i) i 

p01 _(_pO(0)jpla(l) _|_ p0(l)pla(0)^ (gQ) 

(1) i 

T 02 = __L (F 0(0) j p2o ! (l) + p0(l)p2a(0)^ (51) 

(i) i 

7-03 = -(F° (0) F 3a(1) + F° (1) F 3a(0) ), (52) 

Ho 

t 11 = -[-(F^F la ^ + pmpMo) ) + \pmp a m + p^p am) _} Lh npm Fam ^ (53) 

/io 4 4 

^22 = J_[_fp2(0) p2a(l) + p2(l) p2a(0)j + IfpW ^^(0) + p{0) paP(l)j _ \tf2 p(0) pa/3(0)^ ^ 

(i) 1 i 

T 33 _ _L[_(_p3(0)_p3a(l) + ^3(1) ^3a(0)) + If pi 1 ) pa/i(0) + p(0) ^a/3(l)^ (55) 

/^■o 4 
W (1) 1 1 

T 12 = T 21 = __L [(j pl(0) j p2 Q (l) + £1(1)^20,(0)) + ^ 12^0)^/3(0)^ (5g) 

A*o 4 
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(1) (1) 1 

T 13 = T 31 = -[i^ {0) F 3a{1) + Fa {1) F 3a{0) ], (57) 



T 23 = T 32 = ^_^2(0)_p3a(l) + _p2(l) p3a(0)^ ^ 



(1) i (1) 

where T 00 expresses the first-order perturbation to the energy density of the EM fields, T 01 , 

(2) (3) 

T 02 and T 03 indicate the first-order perturbations to the power flux densities of the EM fields 

(i) (i) (i) (i) (i) (i) 
in the x-, y- and z- directions, respectively, while T , T 22 , T 33 , T 12 , T 13 and T 23 represent 

the first-order perturbations to the momentum flux density components of the EM fields. 

By using Eqs. (127j) - (15Dl) and (I4"ll) - (j4"51) . we can calculate the first-order PPFs produced 

(i) 

by the HFRGW. We shall focus our attention to the 01-component T 01 [see Eq. (15 Op ] of the 

first-order perturbation: it expresses the x-component of the power flux density (Poynting 

vector) of the EM fields, i.e., the first-order perturbative power flux density generated by 

the coherent modulation of the preexisting x-component of background power flux. Thus, 

(i) 

the corresponding first-order PPF will be c/hu e T 01 . In this case, although we do not know 
value of the initial phase of "the resonant monochromatic component" of the HFRGW in 
the laboratory frame of reference due to its random distribution, setting the phase difference 
6 = 7r/2 will always be possible by regulating the phase of the GB. The x-component of 
PPF generated by the coherent synchro-resonance (u e = u g ) between the perturbative 
EM fields, Eq. f )4T|) and the GB, Eqs. fl27|) - fl29|) . can then be expressed in the following form 
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n 



(i) 



(i) 

T oi 



/ F 0(0)pia(l) , ^,0(l)_pla(0)\ 



j_,chpx_ _ dipy 
uo e dy dx 



U)e=OJg 



A t ^B { y ) i) k g y(z + Zi 



^ e 4/i [l + (^//) 2 ]V2( 2 + / 2 /; 



sin I — tan 



2R 



fj ' 2^W 2 [1 + (z/f) 2 } 3/2 



cos 



/c 9 r 2 



— tan" 



./ 



exp 



W 2 



W 2 J 4fi R[l + {z/f) 2 ] 1 ' 2 



Fi(y) sin ( ^ - tan" 1 - ) + F 2 (y) cos ' 



2R ~ tan_ 7 



+ 



W 2 



R 2 W 4 



-F 2 {y) sin 



x 



A^^jz + h) 

4/i [l + (z//)2]V2 



Fi(y) cos 



m- 4 ™ 7 



fXigJu 

~2R 



— tan" 



./ 



cxp 



w 2 



(59) 



where 

F 1 (y) = /exp(-^)oos(^)dy, 



(60) 

F 2 ( 2/ ) = /exp(-^)sin(^)^, 



are the quasi-probability integrals. 
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FIG. 1: m > (signal) and ^(background) in 1st (x,y,z>0), 2nd (x<0, y,z>0), 3rd (x,y<0,z>0) 
and 4th (x>0,y<0,z>0) octants. and propagate along opposite directions in the regions 
of 1st and 3rd octants, while they have the same propagating directions in the region of 2nd and 
4th octants. 




FIG. 2: Schematic diagram of strength distribution of and iV^ in 1st and 2nd octants. We 



take note of that 
the yz-plane. 



(0) 



while 



x=0 



x=0 



ivi 1} 



and is "the outgoing wave" to 
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FIG. 3: N£ ' and N£ ' in the 5th (x,y>0,z<0), 6th (x<0, y>0, z<0), 7th (x,y,z<0) and 8th 
(x>0,y,z>0) octants. N x and N x propagate along opposite directions in the regions of 6th and 
8th octants, while they have the same propagating directions in the regions of 5th and 7th octants. 



ATM 1 ' 
4 




FIG. 4: Schematic diagram of strength distribution of A^J and in the 5th and 6th octants. 



Also, we take note of that 



(0) 



x=0 



while 



N 



(i) 



x=0 



N 



(i) 



Unlike Fig. 2, here N x °^ is 



"the imploding wave" to the yz-plane. 



It is very interesting to compare n x °\ Eq. (I31~j) . and rii 1 , Eq. (159|) . From Eqs. (!28|) 
and (l3Tj) . we can see that Ey ^ = at the surface x=0, thus n4 = 0; while numerical 
calculation shows that =0 has maximum. This means that any photon measured 

by a detector (a high-sensitivity microwave receiver) from n x 1 \ x =o will be a signal of the 
EM perturbation produced by the GW. Nevertheless, in the regions of x ^ 0, we have 
n x °^ 7^ 0. At first sight n x will be swamped by the background n x , so that n x has no 
observable effect in this region. However, it will be shown that rii^ and n x °^ propagate 



along the opposite directions in some local regions, and they have the different rates of 
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decay. Thus ni and can be separated by the special fractal membranes (see below), 
so that ni^ (signal), in principle, would be observable. The total PPF passing through a 
certain "typical receiving surface" As at the ?/z-plane will be 

N® = // n^\ x=0 dydz. (61) 



Notice that N x is a unique non- vanishing photon flux passing through the surface i.e., a 
number of photons per second. Eqs. (|59|) and (160]) show that Tlx is an even function of 
the coordinates x, thus n x has the same propagating direction in the regions of x > and 
x < 0; and at the same time, n x is an odd function of the coordinate y, so the propagating 
directions of n x are anti-symmetric to the regions of y > and y < 0(such property 
ensured conservation of the total momentum in the coherent resonance interaction). 
Considering the outgoing (and imploding, i.e., they go in both directions) property of N x ^ 
in the region z > 0(and z < 0) (this is a typical property of the GB [55(), it can be seen 
that and N x ^ propagate along opposite directions in the regions of 1st (x,y,z > 0), 
3rd (x,y < 0,z > 0), 6th (x < 0,y > 0,z < 0) and 8th (x > 0,y,z < 0) octants of 
the reacting region between the magnetic poles, while they have the same propagating 
directions in the regions of 2nd, 4th, 5th and 7th octants, (see FIG1,FIG2,FIG3 and FIG4). 
In our EM system example, all of the following parameters are chosen to exhibit values 
that can be realized in the proposed laboratory experiments that is, they are state of the art: 



(1) P=10W, the power of the GB. In this case, ip rs 1.26 x lO^m" 1 for the GB of the 
spot radius Wo = 0.05m. 

(2) By ^ = 3T, the strength of the background static magnetic field. 

(3) < y < Wo, < z < 0.3m, the integration region As (the receiving surface of the 
PPF) in Eq. ([HI]), i.e., As w 10~ 2 m 2 . 

([!]) I = li + l\ = 0.3m and 6m, the interacting dimensions (or reacting region) between 
the relic GW and the static magnetic field. 

(E]) v e — u g — hGHz(\q = 0.06m), this is typical frequency of the HFRGWs in the 
microwave band Q, J^ v 3 ; Q, 0], and of the HFGW predicted by possible high-energy 
laboratory schemes 
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FIG. 5: The perturbative photon flux (s generated by the HFGW of h 



and v = 5GHz, here detecting bandwidths Av = 1Hz 
x = 0, we take note of that the background photon flux iV : 



(0) 



(i) 



(i) 



ms = W- M /VHz 
8.21 x 10 2 s" 1 at 



x=0 



0[see, Eqs.([28D and (j3T|) ]. thus 



would be an observable value, and and propagate along opposite directions in the 



first octant. 




TABLE I: The x-component of PPFs and relevant parameters. Here A is the root-means square 
value of the HFGW amplitudes, I is interacting dimensions between the HFGWs and the static 
magnetic field, A^and are the x-components of BPF and PPF, respectively. 



A(Hz)~^ 


I = h + him) 
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8.21 x 10 10 


1.24 x 10 22 


3.54 x 10 10 
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8.21 x 10 8 


1.24 x 10 22 


3.54 x 10 8 


1Q -26 
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8.21 x 10 6 


1.24 x 10 22 


3.54 x 10 6 


1Q -28 


6 





8.21 x 10 4 


1.24 x 10 22 


3.54 x 10 4 


10 -30 


6 





8.21 x 10 2 


1.24 x 10 22 


3.54 x 10 2 


1Q -32 
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8.21 


1.24 x 10 22 


3.54 


1Q -34 
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1.24 x 10 22 
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FIG 5 gives result of numerical calculation for Nx .FIG 6 is two dimensional description 
of the numerical calculation for the perturbative photo flux density rix , Eq. (1591) . From 
Eqs.()59D, (J5QJ and ivi 0) and A^we obtained in a 1Hz bandwidth are listed in TABLE U 
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FIG. 7: Envelope of the h rms iy) spectrum for a certain parameter condition. The figure is taken 



from Ref.[40][P. Grishchuk, gr-qc/0504018|. The envelop shows that the r.m.s. values of the 
HFRGW amplitudes in the region of 10 8 -10 10 Hz would be ~ 1CT 30 - 10~ 32 /\/flz,roughly. 



In fact, the expected root-mean-square (rms) value h rms in the GHz band of the di- 
mensionless amplitudes by the different cosmo logical models and parameters are quite 

i y, s y, y, y 



5 if ] . According to optimistic estimation, their or- 
Hz, while conservative estimation may be only 



different m, a, la 14, a, lei, i42i 

ders may be /i rms ~10 -29 — 10 -30 ^ 

/i rms ~10~ 34 — lO -35 ^ VHz. Ref. [40] provides a more average estimation for the r.m.s 
value /w s ~lCr 30 - lCT 32 / 4Wz (see FIG 7). Thus in order to detect the HFRGWs in the 
GHz band, the minimal detectable amplitudes of the detecting systems would be h ~ 1CT 30 
or less at least. Moreover, one often estimates the amplitudes of relic GWs by their energy 
spectra, this is useful because it allows us to quickly evaluate the cosmological importance of 



the generated field in a given frequency interval. However, as pointed out by Grishchuk |39| . 
the primary and more universal concept is the amplitude, not the spectrum density. It is 
the field, not its energy density, which is directly measured by the GW detector. Therefore, 
we listed the PPFs under the different amplitude conditions (h srm ~10~ 22 — 10~ 34 ^ VHz) 
in a 1Hz bandwidth in TABLE [H Of course, possible distribution region of the amplitude 
magnitudes of the HFRGWs may be only /i rms ~10 -30 — 10~ 32 ^ V Hz, there are no so strong 
HFRGWs of /i. rms ~10~ 22 — lCr 28 ^/ VHz, but the estimation of the PPFs can display detect- 
ing ability and sensitivity of the EM system in the different amplitude conditions and in the 
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frequency region. 

TABLE |T] shows that the most interesting region would be "typical receiving surface As" 
at the ?/2:-plane (i.e., plane of x = 0), where 

x= 

(e.g., if A = h rms = 1CT 30 / \jHz and / = 6m, then N c 



while iVi 1} 



has a maximum 

x=0 



= 8.21 x lO 2 ^ 1 ). 

x=0 

We emphasize that for the HFRGW and for the constant amplitude plane HFGW, 
even if they have the same amplitude h rms = 10~ 30 j \/Hz and the frequency v = 5GHz, 
their perturbative effects will be different. For the constant amplitude plane HFGW 
propagating along the symmetrical axis z of the GB, it corresponds to a graviton flux of 
N g = 3.77 x 10 16 s~ 1 at the cross section of the waist of the GB (here the minimum spot 
radius of the GB is equal to 5 cm). Unlike the constant amplitude plane HFGW, due to 
the random property of the HFRGWs, they contain every possible propagating direction, 
thus, as mentioned above (Section II), the propagating directions of the relic gravitons are 
nearing of state of isotropy. In this case, only small fraction of the relic gravitons will pass 
through the cross section of the GB. However, the PPF generated by the resonant coherence 
modulation in our EM system is the first-order perturbation rather than the second-order 
perturbation of usual cavity EM response to HFGWs. Therefore, the strength of PPF is 
proportional to the square root \fN~ g of the graviton flux [i.e., it is proportional to the 
amplitude of the GW, see Eq. (1591) ] and not the graviton flux itself N g (i.e., the amplitude 
squared of the GW). In this case, numerical calculation shows that if the deviation angle 
from the z-axis of the propagating direction of the relic graviton flux is less than 10 degrees, 
then its perturbative effect and that of the graviton flux propagating along the positive 
direction of the z-axis (i.e., best resonant direction, see Section V) are nearly the same. 
Consequently, if all relic gravitons propagating along the deviation angle region (8 < 10°) 
and passing through the cross section of the GB are included, then the relic graviton flux 
at the cross section will be N g « 2.89 x 10 14 s _1 at least. This means that in this case the 
gap between the PPFs produced by the HFRGW and the constant amplitude HFGW will 
be about 1-2 orders of magnitude: this is satisfactory. Notice that then ratio of the square 
roots of such graviton fluxes will be 



N g relic GW = / 2-89 X 10 " _ . ^ , n _ 2 
N g plane GW V 3.77 X 10 



1(1 



8.76 x 10" 2 . (62) 
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(1) 



From Eqs. (|27p - (|29p . (HIT) and (|51|) . we can calculate the 02-component T 02 of the first-order 

(i) 

perturbation, and the corresponding PPF will be c/hu e T 02 , it expresses the first-order PPF 
density n y propagating along the y-direction. By using the similar means, we get n y as 
follows: 



(i) 



a) 

<T 02 > 



Ul e =Ulg 



.__ < P0(0)p2a(l) , p0(l)p2a(0) > 



J_ < LeWb^ > 

\ 4/ i [l + z//]M 2 + / 2 /; 



2/johuj, 



1 Refe)*^* ^ 



9?/ 



dx 



UJ e =Ulg 



sin 



fc s r 2 



tan 



*i(j/)cos( 



tan 



exp 



/ 



cos 



W 2 



k g r 2 
~2R 



+ 



tan 



A®B$%y jz + h) 
fJ 2 f i W 2 [l + (z/f) 2 }^ 



4x 2 \ A®Bi 0) k g (z + h 



Tiuj,, 



1 - 



r 



W 2 J ^ R{l + ( z /f) 



2ii 

2 



- F 2 (y) sin 



^ / >•» 2 

~ tan 7 



(63) 



VT 2 



+ 



/.• 2 
i? 2 



4 

W 1 



x 



4/iq [l + (*//) 2 ]* 



-tan 1 - ) + F 2 (y) cos 



sin I — - — 
1 2R 



— i 

tan — 



exp 



PF 2 



we can see that (1) n y is also an even function of the 



2R f 

By comparing with Eqs. (159]) and (1531 . we can see that (1) rij, 1 '' 
coordinates x and an odd function of the coordinates y. Thus n y has the same propagating 
direction in the regions of x >0 and x<0, and n y 1 has opposite propagating direction in the 



regions of y>0 and y<0. However, unlike property of n x 



(i) 



n. 



(i) 



x=0 



x max j Tly 



(0) 



and we 



also have n 



(i) 



0. Therefore, nip and n y v> have very similar distribution and behavior. 



y=0 



,(°) 



J/=0 



In other words, in almost all regions, will be swamped by the background n K y ' , i.e., n 
has no observable effect. (2) In our case, n£ depends only on the state of ©polarization 
of the HFRGW and it is independent of the state of ©polarization of the HFRGW [see, 
Eq. (|59p ], while riy depends only on the state of ©polarization state and is independents 
of the state of ©polarization [see, Eq. (1631) ]. Thus, the state of polarization displayable in 
the EM system will be only the ©polarization component of the HFRGW rather than the 
©polarization component. 

The quantum picture of the above-mentioned process can be described as the resonant 
interaction of the photons with the gravitons in a background of virtual photons (the statistic 



(0) 
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magnetic field) as a catalyst [h! l^ . i.e., the interaction involving elastic scattering of 
the photons by the gravitons in the background of virtual photons (in the reacting region 
between the magnet poles), which can greatly increase the interaction cross section between 
the photons and the gravitons. In other words, the interaction may effectively change the 
physical behavior (e.g., propagating direction, distribution, polarization, and phase) of the 
partial photons in the local regions, and it does not require the resonant conversion of 
;he gravitons to the photons, the latter corresponds to an extremely small conversion rate 



11] . Consequently, even if the net increase of the photon number (the EM energy) of the 



entire EM system approaches zero, then one still might find an observable effect. In this 
case the requirements of relative parameters can be greatly relaxed, such properties may 
be very useful in order to detect the very weak signal of the HFRGWs. In the case of 
astrophysical phenomenon, an analogous example is deflection of light (an EM wave beam) 
in a gravitational field, which causes the deflection of the propagating direction of the light 
ray, and although there is no any change of the photon number, there is an observable effect. 
Of course in this process the interacting gravitational fields are static (e.g., the gravitational 
field of the Sun). Thus there is no the frequency resonant effect between the GWs and 
the EM waves and the space accumulation effect caused by the coherent interaction of 
the two kinds of waves in the propagating direction, but huge celestial gravitational fields 
compensate for such a shortcoming. In our system the change of the propagating directions 
and distribution of the partial photons in the local regions is caused by the GW, while 
the strong background static magnetic field provides a catalyst to enhance the resonant 
effect between the EM wave (the photon flux) and the GW (gravitons), whose coupling 
compensates in part for the weakness of the HFRGWs. 



V. THE SELECTION OF THE PERTURB ATIVE PHOTON FLUXES 

Because of the random property of the relic GWs, detection of the relic GWs will be more 
difficult than that of the constant amplitude plane GWs. However, we shall show that only 
the relic GW component propagating along the positive direction of the z-axis can generate 
optimal resonant response. It is true that for the relic GW components propagating along 
the x—, y— axes and negative direction of the z-axis, even if u g = u e , the PPFs produced 
by them will be much less than that generated by the relic GW component propagating 
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FIG. 8: The z-axis is the symmetrical axis of the Gaussian beam, k g represents the propagating 
direction of the arbitrary component of the relic GW. 

along the positive direction of the z-axis. Thus the perturbations produced by the relic GW 
components propagating along the different directions cannot be counteracted. In FIG 8 we 
draw the symmetrical axis (the 2-axis) of the Gaussian beam and the propagating directions 
k g of the arbitrary component of the relic GWs. 

In order to compare the PPFs generated by the different components of the HFRGW, 
we shall discuss the perturbations caused by the HFRGW's components propagating along 
some typical directions. 

A. The PPFs generated by the HFRGW components propagating along different 
directions 

Here we assume A = h rms = 1CT 30 ^ y/Hz, v e = v a = 5GHz, the detecting bandwidth is 
one Hz. 

(a) 6 = 0, i.e., the HFRGW component propagates along the positive direction of the 
z-axis. As is calculated: the PPF generated by the component may reach up to 8.21 x I0 2 s~ 1 
in a surface of 10~ 2 m 2 area, (see, TABLE [I]). 

(b) 9 = it, i.e., the HFRGW component propagates along the negative direction of the 



25 



z-axis. 



By using the similar means, one finds 
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(64) 



Different from Eq. (l59i) . each and all terms in Eq. ( 1641) contain oscillating factor 2fc 5 z. We 



TABLE II: The PPFs generated by the resonant HFGW components propagating along the differ- 
ent directions, here = 3T,A®,A & ~ 10~ 30 /VWz, v g = 5GHz, l 2 + h = 6m, 



Propagating directions of the N x (s 1 ) 

resonant components of the 

relic HFGWs 

z 8.21 x 10 2 

-z 2.04x10 

x 4. 07X10- 1 

y 



emphasize that 2kz ~ 209^ for the high-frequency relic GW of u g = 5GHz, the factor 2kz 
will play a major role in the region of the effective coherent resonance. In other words, the 
sign of will oscillate quickly and quasi-periodically change as the coordinate z in the 



26 



region increases. Thus the total effective PPF passing through a certain "typical receiving 
surface" will be much less than that generated by the relic GW component propagating 
along the positive direction of the z-axis, (see Eq. ( 1591) and TABLE UT|) . 

(c) 9 = 7r/2, 0=0, i.e., the propagating direction of the relic GW component is 
not only perpendicular to the symmetrical z-axis of the GB, but also perpendicular to 
the static magnetic field By ^ directed along the y-axis so that it is along the x-axis. 
Here we assume that the dimension of the x-direction of B y ^ is localized in the region 
—^3 < x < I4. Utilizing the similar means the first-order perturbative EM fields gen- 
erated by the direct interaction of the relic GW with the static magnetic field can be given by 



BID 

m 



-A^B^kgcix + l 3 ) exp[i(k g x - u g 



£)] + ^A e B^cexp[i(k g x + u g t)], 



-A^B^kgix + l 3 ) exp[i(k g x - u g t)] - -A^B^ exp[i(k g x + u g t)], 

1 i 
--A^B^k g c(x + l 3 ) exp[i(k g x - u g t)] + -A®Bf>cexp[i(k g x + uj g t)}, 

1 i 

-A®Bfkg(X + l 3 ) eXp[l(kgX - Ugt)] + ^A®Bf) eXp[i(kgX + Ugt)], 

(-L < x < L) 



(65) 



In this case the coherent synchro-resonance (u e = u g ) between the perturbative fields, Eq. 
(1651) . and the GB can be expressed as the following of PPF density, i.e., 
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(66) 



where B® and B^ are the y— and z— components of the magnetic filed of the GB, 
respectively, the angular brackets denote the average over time. Notice that we choose 
the GB of the transverse electric modes, so = 0. By using the same method, we can 
calculate ni , Eq. fl65l) . For example, first term in Eq. fl66l) can be written as 
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It can be shown that calculation for the 2nd and 3rd terms in Eq. (|66|) is quite similar to first 
term, and they have the same orders of magnitude, we shall not repeat it here. Notice that 
unlike ni^ produced by the relic GW component propagating along the positive direction of 
the z-axis, the phase functions in Eq. ( 1671) contain oscillating factor k g (x — z), and because 
it is always possible to choose l 2 + h ^ h + I3, i.e., the dimension of the z-direction of By 
is much larger than its x-direction dimension. Thus, the PPF expressed by Eq. (159]) will be 
much larger than that represented by Eq. fl66|) (see, TABLE [IT]) • 

(d) 9 = 7r/2, 4> = 7r/2, i.e., the relic GW component propagates along the y-axis, which 
is parallel with the static magnetic field By°\ 

According to the Einstein-Maxwell equations of the weak field, then the perturbation of 
the GW to the static magnetic field vanishes [14, 28], i.e., 

4 1} = 0. (68) 

It is very interesting to compare in Eqs. (I59I) . (I64I) . (!67j) and (1681) . as is shown that 
although they all represent the PPFs propagating along the x-axis, their physical behaviors 
are quite different. In the case of6' = = 7r/2, nip — 0, Eq. (|68|) ; when 9 = n and 
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6 = 7r/2, = 0, the PPFs contain the oscillating factors 2k g z and k g (x — z), respectively 
[see Eqs. fl64|) and fl67|) ]. Only under the condition 9 = 0, does the PPF, Eq. (1591) . not contain 
any oscillating factor, but only a slow variation function in the z direction. This means that 
rix produced by the relic GW component propagating along the positive direction of the z- 
axis, has the best space accumulation effect (see TABLE [XT]) • Thus, as previously mentioned, 
our EM system would be very sensitive to the propagating directions of the relic GWs. In 
other words the EM system has a strong selection capability to the resonant components 
from the stochastic relic GW background. Therefore, if real relic GW background has a 
small deviation to the isotropy of space, then it should be possible to provide an HFRGW 
map of the celestial sphere by changing the direction of the symmetrical axis of the GB, or 
alternatively by the utilization of multiple EM detectors. 



B. The separation of the PPFs (signal) from the BPFs 



The reflecting 




The reilecling 
fractal membrane 



FIG. 9: ivi 0) , ivi 1) and N X L) in the 1st and 3rd octants. After N X L) is reflected by the fractal 



(l) 



membrane, (e.g., N^p in the 1st and 3rd octants), N& } and N^ 1 ' will have the same propagating 

(IV 

direction. However, can keep its strength invariant within one meter to the membrane (see, 



(i)' 



r(0) 



e.g., Refs. 



56 



57l|) . while decays as the typical off-axis (radial distance r) Gaussian decay 



rate exp(— 2r 2 /W 2 ) [see, Eq. (j3Tj) ] and attenuated further by superconducting baffles, then the ratio 



Nx /N^ would be larger than one in the whole region of 0.35m < x < lm, although N x ' > N : 
in the region of < x < 0.35m. 



(o) 



In recent years new types of fractal membranes have been successfully developed 



57j, 



581 ] . Firstly, these fractal membranes can provide nearly total reflection for the EM waves 
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Hi': 

The PPFr-i 

detector vj j\fW 



The transmitting 
fractal membrane 




/ 

The transmitting 
fractal membrane 



FIG. 10: ivi 0) , ivi 1} and nP' in the 2nd and 4th octants. Unlike Fig.9, here nP' is the PPF 
transmitted by the transmitting fractal membrane, then the PPF detectors should be put in the 
2nd and 4th octants 



TABLE III: Comparison of the PPF reflected or transmitted by the fractal membrane and the 
BPF in the x-direction, here 5 (0) = 3T, h rms ~ lCT 30 //^:, v g = 5GHz, h + h = 6m and 
detecting bandwidth Au = 1Hz. The PPF reflected or transmitted (defined as )by the 
fractal membrane can nearly keep its strength invariant nearly within one meter distance from 
the membrane 



56 



57, 



58j (or even more attenuated by superconductor baffles along the x-axis). 



Even if according to most conservative estimation to the fractal membranes [68], the photon flux 
reflected or transmitted by the fractal membranes can keep ninety percent of its strength at the 
position of one meter distance from the fractal membranes. Here our comparison is just from this 
conservative estimation. Thus the and would have a comparable order of magnitude in 



the region 35cm< x < 37cm (over a diffraction-limited spot area of ~3xl0 m 2 ). 



The distance to the fractal 
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(photon flux) with certain frequencies in the GHz band; at the same time, they can provide 
nearly total transmission for the photon fluxes with other frequencies in the GHz band 
(the fractal- membrane pattern can be "significantly sub- wavelength in all dimensions" 56) ). 
Secondly, the photon fluxes reflected and transmitted by the fractal membranes can keep 
their strength invariant within the distance of one meter from the fractal-membrane's surface, 
especially if the fractal-membrane reflectors are back-to-back very-shallow (or segmented) 
paraboloid mirrors that focus the PPF on the detectors situated out along opposite ends 
of the x-axis. In this case the diffracted focus spot at each detector exhibits a radius 
of \ g /iT = 6cm/7r ~1.91cm (area of ~ 10~ 4 m 2 ) [9,10]. Thirdly, such frequencies can be 
regulated in the GHz band. Since N x ^ (signal) and N x °^ (background) propagate along the 
negative and positive directions of the x-axis in the first octant (the region of x,y,z > 0), 
respectively, i.e., N^ propagates along the direction toward the fractal membrane, while 
propagates along the direction away from the fractal membrane (see FIG 9). Using the 
reflecting fractal membranes with their plane or paraboloid faces normal to the x-axis, it 
will reflect only A^ 1 and not N^ . Once N^ is reflected (defined as N^ )it will have the 
same propagating direction as . However, after is reflected, it can keep its strength 
invariant within one meter distance from the fractal membrane [56, 57], while decays 



2r 5 

rM'Mr(o) 



as the typical Gaussian decay rate exp(— |^)[see, Eq. (|3T!) ] to each side of the GB (x and y 



directed), then the ratio N x /N x (the signal-to-background noise ratio in the x-direction) 
would be larger than one in the whole region of 0.35m < x < lm (see, TABLE IHIl x is 
the distance from the detectors to the fractal membranes). TABLE ll III shows that the BPF 
Nx°^ is much larger than the PPF Nx in the region < x < 35cm, while the A^and 
have the same order of magnitude at x = 35.09cm, and Nx would be larger than N^ 
in the region of x > 35cm. In other words, in the region the signal-to-background noise 
ratio Nx^ j N x °hn the x-direction might gain up to a comparable order of magnitude. It 
appears better to use the transmitting fractal membranes, because the PPF transmitted by 
the fractal membrane can also keep its strength invariant within one meter to the membrane, 
and the PPF does not change its propagating direction. In this case, the PPF detectors in 
1st and 3rd octants in FIG 9 should be replaced by the detectors in 2nd and 4th octants in 
FIG 10. 

In fact, the circular polarized "monochromatic component" ,Eq.(35), is often called the 
right-handed circular polarization, while a left-handed circular polarized component has 
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following form 



h @ = hu = -h 2 2 = A® exp [i (k g z - u g t)] , 

h® = h 12 = h 2 i = -iA® exp [i (k g z - uj g t)\ , (69) 

where A e ,A® ps A (k g ) /a(t). In our EM system, according to Eqs.(35),(38),(39) and (69), 
the propagating direction of depends on the choice of circular polarization. Thus, if 
the interacting "monochromatic component" is only the left-handed circular polarized state, 
Eq.(69), then the propagating direction of iV"i will be opposite to that generated by the 
right-handed circular polarized component, Eq.(35), and then and propagate along 
opposite directions in the regions of 2nd, 4th, 5th and 7th octants, while they have the same 
propagating direction in the regions of 1st, 3rd, 6th and 8th octants. In such a case, the 
distinguishable PPF from the BPF would be in the regions of 2nd and 4th octants but 
not in the regions of 1st and 3rd octants. 

If the both circular polarizations exist at the same time (In this case the two polarized 
states often have a certain phase difference. More detailed investigation for the issues will 
be done elsewhere), then the PPFs (here we defined them as N^j and N^ z , respectively) 
generated by the right- and left-handed polarized circular components will propagate along 
the opposite directions in the every octant. One of them propagates along the positive 
direction in the x-axis, and another one the negative direction in the x-axis. However,because 
acting effects of fractal membranes to N$ and N^j are quite different, one kind of the two 
PPFs (N^j or N^j) could be distinguished from the BPF. 

For example, in the first-octant (the region of x,y,z>0) N$ and N^j propagate along the 
negative and positive directions in the x-axis, respectively. This means that N^j propagates 
along the direction toward the fractal membrane, while N^u and propagate along the 
direction away from the fractal membrane (see also, Figs. 1, 9, and 10). In this case the 
PPF reflected (or transmitted) by the fractal membrane will be only N$ but not 
and Ni°\ Once is reflected (or transmitted) by the fractal membrane, it will keep 

its strength invariant within one meter distance from the fractal membrane [56, 57], while 
Nx°^ decays as the typical Gaussian decay rate exp(— |^), decay as the exp(— ^) [see 
also, Eq.(59)]. Therefore, the ratio N^/N^ would has a comparable order of magnitude 
in the distance of 35cm<x<37cm from the fractal membrane (see also, Table III). In this 
case, in principle, N^j can still be distinguished from Ni°\ while the other one of the PPFs 
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N^jj will be swamped by N x due to the same propagating direction and the similar decay 
way of them. There is a similar property in the second octant (the region of x<0, y, z>0), 
unique difference is that where N^fj and propagate along the opposite directions, while 
N$> and iVi 0) propagate along the same direction. Thus the distinguishable PPF from the 
BPF would be only N^jj but not N xI . Utilizing the similar means it can be shown that 
in the 1st, 2nd, 3rd and 4th octants, the distinguishable PPF from the BPF will be one 
kind of the PPFs, namely, N x j or . Consequently, role of the fractal membranes looks 
like a "one-way valve" with strong focusing function to the photon flux in the GHz band. 
This property will be very useful to distinction and displaying the PPFs generated by the 
stochastic HFRGW background. 

Of course, if considering other possible noise sources and diffraction effects, the values 
listed in TABLE [J will be further reduced, thus an obvious gap still exists between the 
theoretical schemes and reality. 



VI. THE THERMAL NOISE AND THE EM NOISE. 

At the moment there are no operating prototypes of the EM detecting system, although 
relevant researches and construction of the EM detecting system are already in progress, it 
is difficult to give a complete description for the noise issues. However, since our purpose is 
display and detection of the PPF of about z/=5GHz in the terminal microwave receiver, our 
attention will be focused into two key aspects: 

(1) . What are strength and physical behavior of the PPF (signal) and the BPF (back- 
ground) reaching the microwave receiver; 

(2) . How to distinguish the PPF and other photons caused by noise, such as the thermal 
noise, background noise and external EM noise. Here we shall give a very brief and rough 
review. 

Except for the background photon noise issue just mentioned, there are the thermal 
noise sources and possible external EM noise sources. Because the frequency of the PPF 
(signal) is roughly 5GHz, if the system is cooled down to KT < hu e (K is Boltzmann's 
constant, u e = 27rz/ e , u e = 5GHz), i.e., T < hu e /K ~ 0.24K, then the frequency v m of the 
thermal photons will be less than the v e of the PPF. If the apparatus is kept to a lower 
temperature, e.g., T < 0.024A" or 24mA" (this is well within the current technology), then 
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we have u m ~ 10 2 z^ e . Thus the difference in the frequency band for such two kinds of 
photons would be very great, i.e., the signal photon flux and the thermal photons can be 
easily distinguished. In other words, practically speaking there are no thermal photons at 
5GHz, and in this way the thermal noise can be suppressed as long as the EM detector 
can select the correct frequency. Note that the low temperature is very convenient for the 
operation of the superconductors and the strong static magnetic field. 

For the possible external EM noise sources, using a Faraday cage or shielding covers made 



from such fractal membranes 



56 



571 . |58| | , or from a tight mosaic of superconductor chips 



on the inside surface of the detector's cryogenic containment vessel, would be very effective. 
Moreover, a good "microwave darkroom" can provide an effective shielding environment, 
and in this case possible dielectric dissipation (using a vacuum operation) can be effectively 
suppressed. In this case one would obtain a suitable environment for a measurable signal- 
to-noise ratio. 

Also, the superposition of the relic GW stochastic components will cause the fluctuation 
of the PPFs, even if such "monochromatic components" all satisfy the frequency resonant 
condition (u e = uj g ). However, Eqs. ( I59l . ( 16T1) . ( 164"1) . ( 1671) and ( |68j) show that the metric 
perturbation only influences the strength fluctuation of the PPFs and does not influence 
the "direction resonance." That is, it does not influence the selection capability of the EM 
system to the propagating directions of the relic GWs, and it does not influence average 
effect over time of the PPFs. 

In addition, the values of the PPFs discussed in the present paper depend on the strength 
of the HFRGWs in the GHz band expected by the QIM and other relevant string cosmology 
scenarios (e.g., see Refs, . Because such models and scenarios are somewhat 

controversial, we cannot know in advance how accurate these models and scenarios might 
be. If the strength of the real HFRGWs in the GHz band are much less than the magnitude 
expected by such models and scenarios, even if the required conditions can be satisfied and 
one might still not be able to detect and measure such HFRGWs, then the HFRGW models 
will be corrected. Thus, this scheme might provide an indirect way to test such models and 



scenarios, that is, as suggested by Brustein et al. [631 ]. a null experiment would be valuable. 
In any event, the HFGW generator and detector experiment described in Ref. 9, [loj], which 
operate at about the same frequencies as the HFRGWs, will prove the concept of the present 
detector independently from cosmological experiments. 
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Moreover, there are some issues and problems need further investigation. For example, 
how to generate a typical and high-quality GB, how to suppress distortion of the spot radius 
of the GB and align it, what is concrete correction to the PPF caused the higher order 
modes of the GB, how to further estimate and analysis the relevant noise sources, what 
are concrete influence and correction of the fractal membranes (or plates) to the GB itself, 
how to estimate and effectively suppress diffraction effect by new materials, (e.g. the fractal 
membranes), how to ensure a good vacuum to avoid the scattering of photons and dielectric 
dissipation caused by the dust and other particles, etc. All these issues and problems need 
careful theoretical and experimental study. More detailed investigation concerning such 
issues will be an object of further research and will be studied elsewhere. 

VII. CONCLUDING REMARKS 

1. Although usual analytic expressions of the relic GWs are often complicated, the high- 
frequency asymptotic behavior of them in the microwave band can be expressed as 
simpler forms, and they can be described as superposition of all quasi-monochromatic 
components. The energy density of the HFRGWs is positive definite, and their mo- 
mentum densities have reasonable physical behavior, the EM resonant response of the 
HFRGWs in the laboratory frame of reference can be treated as resonance interaction 
of the quasi-monochromatic HFGWs with the EM fields. 

2. Under the synchroresonance condition, coherent modulation of the HFRGW to the 
preexisting transverse BPFs would produce the transverse PPFs, the PPFs propagat- 
ing along two orthogonal directions of the double transverse polarized electric modes 
of the GB are generated by the pure ® polarization and the pure ©polarization states 
of the HFGW, respectively. The former has maximum at the longitudinal symmetrical 
surface of the GB where the transverse BPF vanishes, but the later and the BPF have 
the same distribution. Thus, the former may provide an observable effect while the 
latter will be swamped by the BPF. 

3. The PPF reflected or transmitted by the fractal membranes exhibits a very small decay 
compared with the much stronger BPF. In our case this is the PPF produced by the 
pure <S> polarization state of the HFGW. Another interesting area would be the region 
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in which the PPF and the BPF might reach up a comparable order of magnitude. 

4. Although an obvious gap still exists between the theoretical estimation and reality, 
there are a potential advanced space and new ways [64-67] to further improve the 
sensitivity and the detecting ability of the EM system. These new ways and technol- 
ogy will include generation of super-strong static magnetic fields (e.g., use of crystal 
channel effect), ultra-high sensitivity microwave single photon detectors such as a cir- 
cuit Quantum Electrodynamics device (CQED) photon detector, Ryberg Atom Cavity 
detector, SQUID array mux, Josephon Junction Arrays, etc., and possible optimized 
combination of them. They are possible to further narrow such gap and provide new 
promise. 
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